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Abstract. We consider the adiabatic demagnetization in the rotating reference frame 
(ADRF) of a system of dipolar coupled nuclear spins s = 1/2 in the external magnetic 
field. The demagnetization starts with the offset of the external magnetic field (in 
frequency units) from the Larmor frequency being several times greater than the local 
dipolar field. For different subsystem sizes, we have found from numerical simulations 
the temperatures at which subsystems of a one-dimensional nine-spin chain and a 
plane nine-spin cluster become entangled. These temperatures are of the order of 
microkelvins and are almost independent of the subsystem size. There is a weak 
dependence of the temperature on the space dimension of the system. 
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1. Introduction 

Entangled states are very important for quantum computing, teleportation, 
cryptography . Although entanglement is a profound concept of quantum information 
theory its experimental realization in many-body systems is an unsolved problem up to 
now. Meanwhile the existing criterions of the existence of entanglement [21 [3] allow us to 
investigate entangled states in conventional NMR experiments [1] . It is well known that 
entanglement emerges in systems of nuclear spins at microkelvin temperatures ^ E]. 
Such temperatures can be achieved with the adiabatic demagnetization in the rotating 
reference frame (ADRF) [7]. 

In the present article we consider the emergence of entanglement in the course 
of the ADRF in the system of nuclear spins coupled by the dipole-dipole interaction 
(DDI). We perform numerical calculations for a chain consisting of nine spins and a 
plane nine-spin cluster in order to prove the emergence of entanglement of its different 
subsystems. We use the Wootters criterion [2] for the investigation of the spin pair 
entanglement. Entanglement of subsystems with larger numbers of spins is investigated 
with the positive partition transposition (PPT) criterion [3]. Entanglement emerges at 
approximately the same temperature of the order of microkelvins for all subsystem sizes. 
However, the temperature depends on the space dimension of the system. 



2. The density matrix of a spin system at the ADRF 

We consider a system of N nuclear spins s = 1/2 coupled with the external magnetic 
field. The system is irradiated by the frequency-modulated radio- frequency (rf) field, 
w{t), which is perpendicular to the permanent magnetic field. The Hamiltonian, Hiab, 
of the system in the laboratory frame can be written as follows 



Hiab = woh + Hdz + 2wi4 cos 



w{t')dt' 







(1) 



where Wq is the Larmor frequency, Wi is the amphtude of the rf field (in frequency units), 
Ina is the projection of the angular momentum operator of the n-th spin (n = 1,2... A^) 
on the a axis (a = x,y,z), la = J2n=i^na, and H^z is the secular part of the DDI 
Hamiltonian [1] which can be written as 

Hdz = ''^Dij{31izljz — I'i Ij), (2) 

i<j 

where Dij is the DDI coupling constant of spins i and j, and /, Ij = lixljx+hyljy+hzljz- 
The spin dynamics of the system is determined by the density matrix, p(t), whose time 
evolution occurs according to the Liouville equation {h = 1) [4j 

^^-^ = [Hu..p{t)]. (3) 

Converting the density matrix, p(t), with the unitary transformation 
pit) = e-''^ Ht')dt'p* ^^yih /* v.{t')de 
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and neglecting the irrelevant terms oscillating with the double Larmor frequency one 
obtains the evolution equation for the density matrix p*{t) 

= [{wo - w{t))h + H,^ + p\t)]. (5) 

According to ([5]) the Hamiltonian, H, of the system can be written in the rotating 
reference frame (RRF) as 

H = A{t)h + Hd, + wih, (6) 

where A(t) = wq — w{t) is the resonance offset of the longitudinal frequency from the 
Larmor frequency. In the course of the ADRF, the offset, A(t), changes slowly in order 
to satisfy to the adiabatic condition [4J 

|A(t)| 



TTwf 



« 1- (7) 



The condition ([7j) means that the resonance offset, A(t), changes so slowly that the 
spin system is in the quasi-equilibrium state at every moment of time [1]. Thus the 
spin system can be described by the thermodynamic quasi-equilibrium density matrix, 

Peq{t), as 

Pe,{t) = e-^^/Z, (8) 

where (3 is proportional to the inverse temperature {/J = h/kT) and Z is the partition 
function. It is important to emphasize that the temperature, T, is the spin temperature 
of the system which is isolated from all other degrees of freedom. The reason of such 
isolation is long spin-lattice relaxation times which exceed spin-spin relaxation times 
by several orders of magnitude [1]. The important consequence for our approach is 
the following. The decoherence effects at entanglement generation with ADRF are 
irrelevant. 

The entropy, S, of the system is given by [1] 

S = -ktT{peg{t)\n[p,g{t)]}. (9) 

We start the ADRF with A > wioc where wioc = {tr {Hj^}/ tr [Z^]}^/^ is the local dipolar 
field. Then the offset, A(t), is the linear function of time 

A = Ao - at (10) 

where Aq and a are the given quantities. Since the entropy, S = const, one can obtain 
the inverse temperature, l3{t), in the course of the ADRF, if the offset, A(t), is known. 
It means that we obtain the density matrix, Peq{t), during the ADRF. In particular, the 
entropy, S, is 

S = kN\n2 + kN\n[cosh ( —]] - -N(3Atanh ( —] (11) 



. 2 , 

at A ^ Wioc- (I lip can be used in order to find the initial inverse temperature of the 
system. 
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3. The reduced density matrix of a spin pair at the ADRF and the 
Wootters criterion 

In order to obtain the reduced density matrix of an arbitrary pair of spins i and j we use 
the approach developed in [5], [6]. The density matrix, Peq(t), of (IE]) can be represented 
as [SIE] 

3 

P= E aitN''xi'®...®x^^, (12) 

where is a number of spins in the system, C,k {k = 1,2,...,N) is one of the 
values {0, 1, 2, 3}, = Ik is the unit matrix of the dimension 2 x 2, = Ikx, 
xl = Iky, xl = Ikz, and ai2^.^Ar^^ is a numerical coefficient. Averaging the density 
matrix of f |T2l) over all spins except spins i and j and taking into account that 
tr{x|*} = {k = 1,2, . . . , N; C,k = 1,2,3) we arrive at the following expression for 
the reduced density matrix, peq\t), of the ith and jth spins 

pifW= E o^'4®^. (13) 

where 

4 = / V • (14) 

The coefficients, of f|T^ can be calculated numerically. Then the reduced density 

matrix of the pair of spins i and j is determined completely. In order to apply the 
Wootters criterion [2j one should find the " spin-fiipped" density matrix, peq\t), which 
is 

pS^it) = {cry ® cry)[p(i^\t)nay ® a,) (15) 
where the asterisk denotes complex conjugation in the standard basis {|00), |01), |10), |11)} 
and the Pauli matrix cr^ = ^ ^ ^ ^ '^^^ calculation of the density matrix, ptq\t), 

and the diagonalization of the matrix product peq\t)peq\t) are performed numerically. 
The concurrence of the two-spin system with the density matrix peq'' (t) is equal to [2] 

C = max{0, 2A — Ai — A2 — A3 — A4}, A = max{\i, A2, A3, A4} (16) 

where Ai, A2, A3, and A4 are the square roots of the eigenvalues of the product 
P^f{t)p^P{t). 

4. The entangled state of a spin subsystem with its environment at the 
ADRF 



The spin pair entangled states are a simple type of entanglement which can be 
described with the Wootters criterion [2] completely. Meanwhile entanglement of bigger 
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subsystems of the system can also emerge in the course of the ADRF. A necessary 
condition of separabihty is described by the Peres criterion [3] which is known as 
the positive partition transposition (PPT) criterion [3]. According to the criterion, 
entanglement emerges, if after the transposition of the density matrix over variables 
of one subsystem there are some negative eigenvalues of the density matrix of the 
system. We wrote a program which allows us to investigate entanglement of an arbitrary 
subsystem and its environment in the nine-spin chain and in the square nine-spin cluster 
in the course of the ADRF on the basis of the criterion p|. The PPT predicts the 
emergence of the spin pair entanglement at the same temperature as the Wootters 
criterion [2]. The PPT criterion uses the sum of the absolute values of negative 
eigenvalues after the partial transposition of the density matrix as the measure of 
entanglement. This measure is called the double negativity [8]. We use the criterions of 
entanglement [21 13] for the investigation of the entangled states in a nine-spin chain and a 
plane nine-spin cluster in the next section. Our choice of nine-spin systems is limited by 
available computational power. Notice that even best available supercomputers do not 
allow one to calculate spin dynamics when the number of spins in the system exceeds 
16 [9]. It is worth to noticing that we consider spin-spin interactions of all spins in 
contrast to the approximation of the nearest neighbor interactions when the number of 
spins in the system can be significantly increased pU?] . 



5. Numerical analysis of entanglement in a nine-spin chain at the ADRF 

The results of the numerical investigation of the spin pair entanglement for a linear 
chain consisting of nine spins coupled by the DDI at the ADRF are represented in 
figure [H We started the ADRF with the offset Aq = lODi where Di is the DDI 




Figure 1. The concurrence, C, versus the dimensionless inverse temperature, f3Di, 
for spins 1 and 2 of the nine-spin chain according to the Wootters criterion. The dash 
line shows the same result obtained with the PPT criterion; wi/Di = 2. 
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coupling constant of nearest neighbors in the chain. The numerical calculations are 
performed for A„ = {10 — n)Di, = 0, 1, . . . , 10. The corresponding values of the inverse 
temperatures, /3„ = h/kTn (n = 0, 1, . . . , 10), are found from ([9]) for the dimensionless 
entropy S/k = 0.5. Figured] shows the concurrence of spins 1 and 2 of the nine-spin chain 
as a function of the dimensionless parameter l3Di. At comparatively high temperatures 
the concurrence is equal to zero (see figure [1]) and the spin system is in a separable 
state. When the temperature is getting sufficiently low in the course of the ADRF 
the concurrence is sharply increasing and entanglement emerges. The entangled state 
appears at PDi ^ 1.1 (see figure [1] ), i. e. at T 0.5fiK when Di = 27110"^ s^^. 
Notice that ordered states of nuclear spins were observed in a CaF2 single crystal at 
microkelvin temperatures [?]. 

0,5 I I 




0,2 0,4 0,G 0,8 1,0 1,2 

Figure 2. The double negativity versus the dimensionless inverse temperature, (iDi, 
for spin 1 of the nine-spin chain and the other spins; wi/Di = 2. 

The dash line in figure [1] demonstrates that the result for the spin pair entanglement 
obtained with the PPT criterion coincide practically with the ones obtained with the 
Wootters criterion. Figure [2] shows the double negativity versus the dimensionless 
parameter j3Di for the first spin of the chain (the first subsystem) and the other spins 
of the chain (the second subsystem) in the course of the ADRF when the dimensionless 
entropy S/k = 0.5. One can conclude that entanglement emerges at (3Di 1.1. This 
result is close to the one for the spin pair entanglement. Figure [3] shows the double 
negativity versus the parameter fiDi for the first three spins of the chain (the first 
subsystem) and the other six spins of the chain (the second subsystem) at the same 
conditions. Here entanglement is getting sufficiently large at (3Di 1.1. In fact we have 
found that entanglement of different subsystem emerges approximately at the same 
temperature. 
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Figure 3. The double negativity versus the dimensionless inverse temperature, f3Di, 
for the first three spins of the nine-spin chain and the other spins; wi/Di ~ 2. 
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Figure 4. The square cluster of nine spins. 



6. Entanglement in the square cluster of nine spins 

The suggested method allows us to study entanglement in the systems of arbitrary space 
dimensions. As an example, we consider the square cluster consisting of nine spins (see 
figure H]). The dipolar coupling constant of spins j and k (the numbers of spins are 
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pointed in figure Hj) is 

D,fc = ^(l-3cos2^,fc), (17) 

where 7 is gyromagnetic ratio, Vjk is the distance between spins j, k and 6jk is the angle 
between the vector, fjk, and the external magnetic field, Hq. The simple analysis yields 

2n ^ 9({(j-l)/3}-{(A:-l)/3})^ 

ilU - l)/3] - [{k - + 9({(j - l)/3} - {(A: - l)/3})2' ^ ^ 

and 

Tjk = a^ilij - l)/3] - [{k - l)/3]y + 9({(j - l)/3} - {(A: - l)/3})2, (19) 

where a is the distance between the nearest neighbors in the cluster, [q] is the integer 
part of q and {q} is the fractional part of q. 




Figure 5. The double negativity versus the dimensionless inverse temperature, /3Di, 
for spins 1, 2, 4 (first subsystem) and spins 5, 6, 8, 9 (second subsystem), in the square 
cluster; wi/a — 2. 

Figure [5] and figure [6] show that the entangled states emerge at PDi = 1.3 in the 
two-dimensional cluster. Although the temperatures of the appearance of entanglement 
are below here than in the one-dimensional case, they are almost the same for different 
subsystems. The processes of destructive interference are more effective in the two- 
dimensional cluster than in the one- dimensional chain. They lead to a loss of many-spin 
correlations which are responsible for the appearance of the entangled states. Thus the 
temperature of the emergence of entanglement is lower in the two-dimensional cluster 
than in the one- dimensional chain. 
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Figure 6. The double negativity versus the dimensionless inverse temperature, /3Di, 
for spins 1, 2 (first subsystem) and spins 3, 6, 9 (second subsystem), at Wi/a = 2 in 
the square cluster of nine spins. 



7. Conclusion 

We investigated numerically entanglement in the chain of nuclear spins and in the square 
cluster with the DDI in the course of the ADRF using a special computer program. We 
showed that the entangled states emerge at microkelvin temperatures for typical DDI 
coupling constants. Two criterions [2|, [3] of entanglement yield the same results for the 
spin pair entangled states. Entanglement of different subsystems emerges approximately 
at the same temperature and the pairwise entanglement can be used as an indicator of 
entanglement of bigger subsystems. It is also worth to notice that we take into account 
the DDI of the remote spins in contrast to the works [5], E] where the nearest-neighbor 
interactions were only considered. The performed calculations show that there are no 
the entangled states of remote spins both in one-dimensional and two-dimensional cases. 
Entanglement emerges only if subsystems are in a direct contact. Entanglement of 
different subsystems at microkelvin temperatures suggest possible applications of linear 
spin chains in quantum information processing. 
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